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Abstract —Modern society critically depends on the services 
electric power provides. Power systems rely on a network of 
power lines and transformers to deliver power from sources of 
power (generators) to the consumers (loads). However, when 
power lines fail (for example, through lightning or natural 
disasters) or when the system is heavily used, the network is 
often unable to fulfill all of the demand for power. While systems 
are vulnerable to these failures. Increasingly, sophisticated 
control devices are being deployed to improve the efficiency 
of power systems. Such devices can also be used to improve the 
resiliency of power systems to failures. In this paper, we focus 
on using FACTS devices in this context. A FACTS device allows 
power grid operators to adjust the Impedance parameters of 
power lines, thereby redistributing flow in the network and 
potentially increasing the amount of power that is supplied. 
Here we develop new approaches for determining the optimal 
parameter settings for FACTS devices in order to supply the 
maximal amount of power when networks are stressed, e.g. 
power line failures and heavy utilization. 

1. Introduction 

One of the challenges facing modern infrastructures is the 
possibility of events that cause large-scale damage. These 
events include hurricanes (i.e. Superstorm Sandy, Hurricane 
Katrina, etc.), ice storms, earthquakes, etc. During these 
events, significant portions of a society may be without 
services, such as electric power, that these systems provide. 
As a result, there is considerable interest in developing 
approaches and methods for mitigating the effects of these 
events and restoring these systems as quickly as possible 
[1], [2]. One approach for mitigating the effects utilizes 
controls in the system to limit the impact of large-scale 
events. While installing such controls solely for the purpose 
of responding to big events may be too expensive in general, 
many systems have experienced a rapid rise in deployment 
of advanced control technologies to improve efficiency, cost, 
and reliability to small-scale events [3]. Here, we focus 
on how existing adopted technologies may be used during 
emergency situations to mitigate the impacts of large-scale 
events. In particular, we focus on how Flexible AC Transmis¬ 
sion System (FACTS) [4] improve power throughput when 
a system is stressed due to damage or increased utilization 
of the system. 

The flow of electric power in power systems is governed 
by complicated non-linear physics. From a practical perspec¬ 
tive, this means power flows along paths of least resistance 
from sources of power to consumers of power. As a result, 
operators have limited ability to direct how power flows in a 
network. This lack of control has a number of consequences, 
but the most important consequence here is that a network’s 
capacity may be underutilized. One solution to this problem 
is FACTS devices. Using a wide variety of technologies. 


FACTS devices allow an operator to modify the resistivity 
of power lines, to shift power away or towards portions of 
a network. As a result, power systems may use more of 
their capacity, are more evenly utilized, and less expensive 
power generation may be dispatched. While their primary 
use is in daily operations to improve efficiency, security, and 
economics, FACTS devices, as shown in this paper, are also 
useful for improving system response during stressed and 
adverse conditions. 

In general, most recent work of FACTS device related 
optimization and control has focused on how FACTS can 
contribute to the stability [5], [6], [7], [8], [9] and the 
maximum load-ability of a network [10], [11], [12], [13]. 
There is also work identifying where to place FACTS devices 
such that the load-ability of a network is maximized [14], 
[15], [16]. [17] provides a survey of different goals and 
methods that are used in FACTS optimization. All of those 
approaches use heuristics (for example genetic algorithms, 
particle swarm) and run simulations on small power networks 
(up to 59 buses). An exception to these papers is the work 
[18] that studies the FACTS placement problem on the 2736 
bus polish network using the Linear DC model. In this paper 
we also use the Linear DC model and in contrast to existing 
work, we study the maximum throughput of a network on the 
big 2736 buses polish network. We also focus on developing 
a globally optimal approach as opposed to heuristics. 

The key contributions of this paper are as follows: 

• An optimization model for maximizing power through¬ 
put in power systems when FACTS devices are avail¬ 
able. 

• Empirical results demonstrating how FACTS devices 
increase power throughput in stressed and damaged 
systems. 

• A discussion of the computational complexity of maxi¬ 
mizing power throughput in power systems with FACTS 
devices. 

• Unlike similar problems, such as transmission switching 
[19], our empirical results suggest that, in practice, the 
problem is often tractable to solve to small optimality 
gaps, even for large systems. 

• An algorithm for improving the computational perfor¬ 
mance of commercial solvers on this problem. 

The remainder of this paper is organized as follows. Section 
II defines the problem. Section III discusses the complexity 
of the problem and Section IV discusses our algorithm for 
solving the problem. Section V describes the results of our 
computational experiments and we conclude with Section VI. 



II. Problem Definition 


Power flows in electrical power networks are described by 
nonlinear, steady-state electrical power flow equations (Alter¬ 
nating Current Model, AC). In this paper we use the Linear 
DC (LDC) model which is an approximation (linearization) 
of the AC model [20]. This model ignores reactive power 
and resistance and assumes that all voltages magnitudes are 
one in the per-unit system. What remains are susceptances 
(the negative inverse of the reactance), the capacities and the 
phase angles of the voltages. The flow on a line is similar 
to that of DC currents. The susceptance is the counterpart 
to the DC resistance, phase angles are the counterpart to the 
DC voltages and the power is the counterpart to the current. 
FACTS devices are physical devices allowing the (otherwise 
constant) susceptance parameter to vary. 

Definition 1: A FACTS Linear DC network (FLDC net¬ 
work) is a tuple M = {N , N q, NE) where N is the set 
of buses', Nq C N is the set of generators', N^ C 
is the set of loads', and E C V 2 {N) x I(K+) x IR+ is the 
set of lines with their susceptance limits and 
A Linear DC network (LDC network) network is an FLDC 
network without FACTS devices, i.e. all susceptances are 
fixed. We define functions B* : E ^ I(]R+) for the suscep¬ 
tance limits and C : E ^ IR+ for the capacities of power 
lines. For a line from a to b with susceptance limit [s, t] and 
capacity p we use notation If the susceptance is fixed, 

i.e., f = s, we write ab^. We may also ignore these values 
and simply refer to the line by ab. While this model does 
not explicitly give upper bounds on the generation or load of 
a bus, such constraints can be modeled by connecting these 
buses to the network through a single line whose capacity is 
the maximum output/intake of the bus. 

We now introduce the notations and equations describing 
FLDC network power flows. We assume a fixed FLDC 
network Af = (TV, Nq,N j^,E). The generation and load at 
a bus are given by functions G : TV —)■ K’*' and L N ^ IR+ 
such that Va € N\Nq : := 0 and Va € N\N: L^ := 

0. Also, we define functions B : E ^ K+ and 0 : TV —)■ K 
such that is the susceptance of line e and 0^ is the 

phase angle at bus a. The flow on a line is given by 
function F : E ^ M.. While the lines of FLDC networks are 
undirected, orientation is needed to describe flows. However, 
the concrete orientation we choose does not influence the 
theory. To that end, whenever we define a line, we abuse the 
notation ab to indicate that > 0 whenever the flow goes 
from a to b and < 0 otherwise. 

The LDC network model imposes two laws: Kirchhoff’s 
conservation law and the LDC network power law. 

Definition 2: A triple {F,G,L) satisfies Kirchhoff’s con¬ 
servation law if Va G W : J^abGE Kb “ Y^baaE Ka = 
Ga - La- 

Definition 3: A triple {B,Q,F) satisfies the LDC net- 



(a) An LDC network. (b) The MPF. 



(c) An FLDC network (d) The MFF. 

network. 


Fig. 1: Examples for MFF and MPF. 


work power law if Va6 G E : = F^t,(0f, — 0^) 

Definition 4: We call a tuple {B,Q,F,G,L) di feasible 
solution if: {F,G,L) satisfies Kirchhoff’s conservation law; 
(F,0,F) satisfies the LDC network power law; Ve G F : 
Ba G B: and Ve G F : |FJ < G,. 

Finally, we use Iff for the set of all feasible solutions of Af. 

The maximum flow AviE{N) of a network A/" is a triple 
{F,G,L) that maximizes the generation w.r.t. respecting 
Kirchhoff’s conservation law; the line capacities and the 
generation and load bounds. We now define two variants of 
this problem for power networks. The maximum FACTS flow 
(MFF) additionally has to satisfy the LDC power law. The 
maximum potential flow (MPF) is a variant of MFF where 
susceptance is fixed (i.e., it applies to an LDC network). 

Definition 5: The maximum FACTS flow (MFF) of 
an FLDC network N is defined as A4EE{N) := 

^ti.X(B^Q^F,G,L)GTjv J2gGN Gg- 

The maximum potential flow (MPF) of an LDC network A/^ 


is defined as AAVF{N) '.= ,gY. geN Gg- 

Figure [T^ shows an LDC network where p is a generator 
(box), I is a load (house) and & is a bus (sphere). We omit 
the susceptance and capacity of a line when its value is 1. 
Here, the MF for this network is 34 whereas in the LDC 
model, we only can supply 16 as shown in Figure [Tb|because 
the congestion of the edge bl\ constrains the phase angle 


(written as A = in the buses) between g and 1. Figure Ic 


shows a variant of the network with two FACTS devices. 
These devices allow the maximum generation to reach 28 as 
shown in Figure 


HI. Computational Complexity 

Finding the solution to the MPF is known to be polynomial 
as it can be described as a linear program (LP). In our 
workshop paper [21] we prove that finding the solution to 
the MFF is NP-complete even for simple network structures. 
For completeness, we discuss the main idea of the proof 
here. We also present three special cases where the problem 
is polynomial: the network is a tree; all lines afophave 


'W.l.o.g. buses with load and generation can be split into separate buses. 
^W.l.o.g. we assume that no two lines connect the same pair of buses. 


^Normally, the susceptance is a negative value and the flow equation is 
— ©(,)■ For notation simplicity, we make the susceptance a positive 
value and multiply the flow equation by —1. 















C=x 


FACTS devices with either s = 0 or i = oo. In all three 
cases the MFF is equal to the MF. 

Lemma 3.1: Let JV be an FLDC network with a tree 
structure then MFF{N) = MT(Af). 

Proof: This is a consequence of the absence of cycles. 
Hence there are no cyclic dependencies on the phase angles 
which allows us to chose them in a way to match any optimal 
solution of the traditional max flow. ■ 

Lemma 3.2: Let JV be an FLDC network where for all 
lines a&j,*’*': s = 0 then J\AFF[N) = J^T{JV). 

Proof: Let {F,G,L) be an acyclic optimal solution 
of the MF. We have to define susceptance B £ B* and 
phase angles 0 such that the DC power law is satisfied. 
First, we define preliminary phase angles 0' using arbitrary 
positive values with the restriction that they respect the flow 
directions, so Vo 6 G i? : 0(j > 0^ if < 0; 0„ < 
if > 0; and 0^^ = 0^, if F^^ = 0. This defines 
susceptances B' via the DC power law: Vab G F : B , := 
, , . We now scale these susceptances such that they fit 

into their limits. Let x := | G F}. By 

^ab ^ 

setting B xB' and 0 := ^0^ we obtain Va 6 G E : 
Fab = 

and hence (5,0, F) satisfies the LDC power law. Also, 
using the definition of x, for an arbitrary G F we 

have X < and hence 0 < < s. ■ 

Lemma 3.^: Let JV be an FLDC network where for all 
lines we have t = oo then J\A.FF{N) — J^JF{N). 

Proof: Similar to the proof of Lemma |3.2| ■ 

In [21] we prove that deciding whether or not 
J^JFJ-{N) > X for some given a; G K is strongly NP- 
complete in general and NP-complete when we restrict the 
network structure to cacti. A network is a cactus if every line 
is part of at most one cycle. The key element of these proofs 
are choice networks. While the workshop paper presents 
the pure mathematical proof of the properties of choice 
networks, the reminder of this section presents the underlying 
idea that makes these properties possible. 

Choice networks are used to encode (discrete) choices 
that characterize NP-hard problems. Given an a; G K"*", a 
choice network can be regarded as a black box with port 
p where we have external generation of either a; or 0 at p 
in order to achieve the inner maximum generation. External 
generation indicates that the black box acts as generator for 
the network that is connected to to p. Inner generation is the 
generation produced inside the black box. Figure [^presents 
our generation-FACTS-choice network GFCN^ 

In the following we describe the way the choice network 
works for the case x = 1. In Figure the generator can 
deliver 1 unit of power. Let the flow to p be re and the flow to 
the network A be r, so that w+r = 1. For the port p we have 
JAFF (p) = w. The MFF of port p depends on w according 
to a linear function of the form f{z) = n + mz. Here, n is 
the base generation and m is the generation ratio. The port p 
has a base generation of 0 and a generation ratio of 1. In the 
DC model and in a network without FACTS device, n and m 
are independent from the input. The usage of FACTS devices 
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Fig. 2: A choice network with port p. 



Fig. 3: Example encoding for {M,S) = ({a, 6 , c, d, e,/}, 
{{a, 6 , c},{ 6 , c,d},{d,e,/}}). 


allows us to construct networks where both values change 
if the input exceeds some threshold. As we will show later, 
network A is constructed in a way that JAFF{A) = 5.1 + |r 
if r < 0.65 and JAFF{A) — (5.1 ~ |) + |?' otherwise. 
Using the equation w = 1 — r and assuming r < 0.65 we 
have JAFF{GFCNp i) = 1 —r+5.1 + |r = 6.1—which 
has a single maximum of 6.1 at r = 0. On the other hand, if 
r > 0.65 then we have JAFF{GFGNp = 1 — r + (5.1 — 
|) + |r = 6.1 — I + I?” which has a single maximum of 6.1 
at r = 1. This implies that there are exactly two solutions 
to achieve a flow of 6.1 and that port p has 0 or 1 units 
of power. Looking at the general case, we see that p either 
generates 0 or x. 

The MEE of A depends on the value of r: JAFF{A) = 
5.1 + |r if r < 0.65 and MFF{A) = (5.1 - 5 ) + |r 
otherwise. The key feature to make it possible that there are 
exactly two solutions for the MEE is that the generation ratio 
changes from a value less than one (|) to a value greater than 
one (|). We now explain why A has two generation ratios. 
The network A consists of two networks in sequence: B and 
C. The network C has a fixed generation ratio of 2. Every 
unit of power that enters d has to go to /. This increases 
the phase angle difference between d and / by one. Hence 
the phase angle difference between e and / increases by 
one and gives us an additional unit of power. Network B 
has two generation ratios: ^ for an input less then 0.65 and 
I otherwise. Because B and C are in sequence, the ratios 
multiply and we get ratios of | and | for A. 

Using choice networks, we reduced the exact cover by 3- 
set problem into the MEE problem [21]. Eigure [re presents 
the idea of the proof. 

Theorem 3.4: Let x G K and JV be an ELDC network. 
Deciding if JAFF{JV) > x is strongly NP-complete. 

Proof: We prove this by reduction from the exact cover 
by 3-set problem. Given a set M and a set of subsets S C 









































V{M) where every element of S has exactly 3 elements, 
decide if there exist a set T C S' such that Uxgt X = M 
and VXi, X 2 G r : Xi 7 ^ X 2 ^ n X 2 = 0. 

For an instance (M, S), we define the FLDC network 
{N,Nc;,N^,E) with Nq := {(?}, 

{1}, N := iVg U U UxGs{^^^} U Ua;GM{4 and 
E •= { 9 ^ 3 } U UxgS Ua:Gx{^^a;i}) U [JxeM{9^lJ^^2}- 
We then define + Yl,x&s GFCN^ and 

we have: MFF{N^’^) = 3 + 18.3|S| + \M\ ^ 
(M, S) is solvable. An example encoding for the exact cover 
problem (M, S) = {{a,b,c,d,ej},{{a,b,c},{b,c,d}, 

{d, e, /}}) can be found in Figure]^ The choice networks are 
represented by 3D boxes where the number 3 in parenthesis 
is the value chosen for x. 

For each b G M we have a bus b. The network is 
constructed such that = 3 + 18.3|S| + \M\ 

iff all choice networks have their inner maximum generation 
(3 • 6.1) and all lines from the generator are congested. 
This implies a phase angle difference of 3 between the 
generator and the load and a phase angle difference of 1 
between b and the generator. Hence we have a phase angle 
difference of 2 between the load and b. For b to satisfy 
Kirchhoff’s conservation law, b receives one unit of flow 
from a GFCN^ with b G X. Because GFCN^ can 
only generate 3 if it wants to generate anything, all c G X 
get one. Hence, we can only achieve the proposed MFF value 
iff (M, S) is solvable. ■ 

IV. Global Optimization 

In this section we present a mixed integer program (MIP) 
model to find the solution to the MFF. The main problem in 
solving the MFF is the right-hand-side of the DC power law 
= ^a 6 (®b ~ ®a)' This Constraint contains the product 
of two variables which, in general, is a non-convex quadratic 
constraint (as implied in the previous section). 

However, for solving the MFF, we do not need to know the 
susceptance values, we only need to ensure that there exists 
valid susceptance values such that the phase angle difference 
and the line flow are bound together. The solution region 
for just the phase angle difference and the flow consists 
of two convex regions. There is one region for positive 
phase angle differences and one region for negative phase 
angle differences. Given an edge 06 ^’*^ the positive region is 
described by s(0j —< F^^ < f(0j —0^). To formulate 
the problem as a MIP, we introduce a binary variable dab 
that represents the choice of phase angle difference direction. 
Assuming that M is some upper bound for the phase angle 
difference, the MIP model can be found in Model [T] 

In this model. Equation describes the objective function 
that maximize the amount of flow in the system by maxi¬ 
mizing the total generation. Equation constrains the flow 
direction variables as binary. Then, equations and force 
either the positive or negative phase angle difference of ab 
to be 0 depending on the choice of d. Equation states that 
the difference in phase angles (right hand side), as stated by 
0 , is equal to the sum of the negative and positive phase 
angle differences as stated by A. This constraint couples the 


maximize 

aeN 

(1) 

subject to 

dab G {0,1}, ab G E 

(2) 


0 < A+ < dabM, abGE 

(3) 


0 < A“^ < (1 - dab)M, ab G E 

(4) 


^tb - Kb = ©6 - ©a. abGE 

(5) 


sA+ < F+ < fA+, a6jr-‘l G E 

(6) 


sA- < F- < fA-, G F 

(7) 


Fab = Ftb - Fab> e F 

(8) 


Fab - E K = G'a - L^, 

ab£E ba£E 

aGN (9) 


— Cab ^ Fab ^ Cab, ab G E 

(10) 


Model 1: MIP model to solve the MEE. 


flow direction (edge based variables) with the phase angles 
(node based variables). Equations and state that the flow 
in the positive or negative direction can be any value that is 
consistent with the phase angle differences and the bounds 
on the susceptances. Equation states the flow on ab is the 
sum of the negative and positive flows. Equation states 
that flow must be balanced at each node. Finally, equation [T0| 
constrains the flow on ab to be smaller than its capacity. After 
this problem is solved, the susceptance of ab is calculated as 

^ab — e^-e„' 

During our experimental testing, we observed cases where 
commerical solvers such as Gurobi and Cplex were unable 
to find provably good solutions on large problems. After 
one CPU hour the optimality gap was larger than 400%, 
in some cases. Based on these observations, we developed 
an approach to boost the performance of the solvers by 
providing high quality initial solutions. We refer to this 
approach as the iterative method f/M).The IM is based on 
two observations. First, finding the MPF is polynomial. That 
means that finding the maximum flow in a network with fixed 
susceptances is easy. Second, if we fix the binary variables 
in the MIP model, we are left with an LP. Fixing the binary 
variables is equivalent to fixing the phase angle difference 
for all lines. We refer to the maximum possible flow possible 
when the sign of the phase angle differences fixed as the 
maximum variable flow (MVF). The IM works as follows, 
given some random valid susceptances {B) we find the MPF. 
Then we take the sign of the phase angle differences {D) 
from this solution, fix them and solve the MVF problem. 
From this solution we take the susceptances {B) and put 
them into the MPF model to start the next iteration. This 
process is described in Algorithm We continue iterating 
until we converge to a solution. 

V. Experimental Results 

For our experimental results we focused on evaluating the 
computation required for solving the MPF, both with and 
without the IM. We also compared the solutions obtained 



Data: susceptances B 

Result: phase angle difference directions D 
MPF = -1; 

MVF = 0; 

while MPF < MVF do 

MPF, D = solveMPF{B)-, 

MVF, B = solveMVFiD); 

end 

Algorithm 1: The iterative method (IM). 
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Fig. 4: IEEE 2736 damage scenarios that remove randomly 
selected lines (x axis). The y-axis shows the maximum 
throughput for MPE and MEE. 


by the MPE with the solutions of MEF in order to provide 
some evidence of the types of solution improvement one 
might expect when FACTS devices are used. We used the 
IEEE test problems provided with Matpower [22] to test 
our approach. In general, on smaller networks, the MEE is 
computationally easy to solve (under 1 CPU minutes) and 
the differences between MPE and MEE are small. Thus, 
we focus on presenting results produced on the IEEE 2736 
bus system (based on Poland’s power grid) due to its large 
size and computational challenges. The results were obtained 
on a computer with an Intel Core i7-4702HQ Quad Core 
processor using Gurobi 5.5 [23]. 

We used two approaches to create variants of the IEEE 
problems that mimic damage and heavy utilization. Our 
first method removes random lines to simulate line failures. 
Eigure shows results on IEEE 2736 for 60 scenarios where 
random sets of lines are removed from the power network. 
The graphs compare the quality of the solution for MPE with 
MEF. MFF was initialized with IM and was allowed to run 
for 10 CPU minutes. In each case we selected a random 
set of lines that have FACTS devices according to a uniform 
distribution. Each device was allowed to vary the suceptances 
30%. With these results, we see that as the amount of damage 
increases, the benefits provided by MEE increase (up to 5%). 

Our second method scales the upper bound of load and 
generation to model increased utilization. We refer to the 
scaling factors as generation and load congestion factors, 
respectively. Table [Ijpresents the solutions for MPE, MEF and 


their improvement percentage for different congestion factors 
when 30% of the lines (randomly chosen) have FACTS 
devices that are allowed to vary from their initial susceptance 
by 30%. MFF was initialized with IM and was allowed to run 
for 10 CPU minutes. The first row and column of this table 
note the generation and load congestion factors. In this table, 
the results indicate that as the congestion factors increase, the 
benefit of MPF over MFF increases by up to 3%. 

In most of the results discussed in this section, Gurobi is 
able to find solutions with a gap of 0.1% in under an hour of 
CPU time without the need for IM. As part of our study, we 
generated a large number of random cases for different levels 
of FACTS penetration and different capabilities of FACTS 
to vary susceptances. In some of these cases, even with 1 
CPU hour available, Gurobi was unable to find reasonable 
solutions. In general, we observed this behavior when the gap 
between MPF and MFF was large (between 6 and 11%). 
We show some results for one case when the generation 
congestion factor was 2.375 and the load congestion factor 
was 2.75 in|n] The MPF of this network is 419.9. The first 
two columns of this table show the percentage of lines with 
FACTS devices and the allowed interval of susceptance (as 
a percentage of the original suspectance). The MFF column 
shows the results of MFF without IM. The columns labeled 
Obj, RT, and MIPGap present the solution, run time in CPU 
seconds, and optimality gap, respectively. The column 3- 
IM presents the results of the iterative method in terms of 
solution quality (IM), run time in CPU seconds (RT), and 
iterations (Calls). Here, IM is executed 3 times with different 
starting solutions (susceptance upper bound, sucsceptance 
lower bound, and suscepatance mid point), and the best result 
is reported. As seen in these results, the IM method produces 
a high quality solution in a short amount of time. In the 3-IM 
+ MFF column, we show the results after 10 CPU minutes 
when IM is used as an initial solution for MFF. Though the 
MIP solver is unable to improve the solutions much, it is 
able to show that the IM method (in most casess) produces 
solutions with substantially small optimality gaps and that 
are significantly better than MPF. 

VI. Conclusion 

In this paper we discussed an optimization model for 
maximizing the throughput of a network utilizing FACTS 
devices and an algorithm that improves the computational 
performance of commercial solvers on this problem. In 
general, even for large scale problems, our approach is able 
to solve the problem to near optimality in under 10 CPU 
minutes. On large scale problems, using FACTS devices 
when the network is over utilized or damaged is able to 
improve the throughput between 2 and 5% (on occasion, 
higher). While this appears to be a modest improvement, 
when discussing damage scenarios that result in millions 
of people without power (such as Hurricane Sandy), a 
5% improvement corresponds to a significant number of 
people. We also provided a discussion that illustrates that 
the problem is NP-complete in general. 



































































































































TABLE I: A comparison of MFF and MPF on IEEE2736 for different congestion factors. 30% of the lines have FACTS 
devices 


Generation 

1.5 

2.0 

2.5 

3.0 1 

Load 

MPF 

MFF 

Improv. 

MPF 

MFF 

Improv. 

MPF 

MFF 

Improv. 

MPF 

MFF 

Improv. 

1.5 

270.56 

271.0 

0.16 

271.12 

271.12 

0.0 

271.12 

271.12 

0.0 

271.12 

271.12 

0.0 

2.0 

303.05 

303.13 

0.03 

347.8 

353.47 

1.63 

354.99 

358.68 

1.04 

358.17 

360.33 

0.6 

3.0 

303.13 

303.13 

0.0 

392.83 

401.49 

2.2 

433.85 

445.88 

2.77 

450.15 

462.18 

2.67 

4.0 

303.13 

303.13 

0.0 

401.26 

403.12 

0.46 

461.17 

473.41 

2.66 

490.41 

502.01 

2.36 


TABLE II; Results comparing MFF with IM on problem IEEE2736 with different FACTS devices. The MPF of the network 
is 419.19. IM significantly improves the quality of solutions. 


Net 

MFF (Ih time limit) 

3-IM 

3-IM + MFF (lOmin time limit) 

Improv. 

FACTS% 

Interval% 

Obj 

RT 

MIPgap 

Obj 

RT 

Calls 

Obj 

RT 

MIPgap 

to MPF in % 

30 

40 

405.32 

3600 

15.69 

463.83 

306 

83.00 

463.86 

600 

1.29 

10.66 

60 

70 

455.22 

3600 

3.07 

466.87 

320 

91.00 

466.98 

600 

0.71 

11.4 

60 

10 

91.47 

3600 

398.14 

445.78 

191 

45.00 

445.80 

600 

2.34 

6.35 

10 

40 

449.72 

3602 

4.25 

458.21 

188 

23.00 

460.21 

600 

2.13 

9.79 

30 

10 

88.85 

3600 

424.37 

457.66 

224 

41.00 

457.69 

600 

1.80 

9.18 

60 

40 

343.99 

3603 

35.85 

460.12 

218 

41.00 

460.15 

600 

1.62 

9.77 

90 

40 

376.41 

3600 

23.58 

451.43 

176 

21.00 

454.08 

600 

2.52 

8.32 

90 

70 

455.22 

3600 

3.07 

466.87 

320 

91.00 

466.98 

600 

0.71 

11.4 

10 

10 

385.63 

3600 

21.42 

463.98 

304 

73.00 

464.04 

600 

1.07 

10.7 

10 

70 

455.22 

3600 

3.07 

466.65 

230 

67.00 

466.73 

600 

0.55 

11.34 

30 

70 

455.22 

3600 

3.07 

466.87 

254 

65.00 

466.98 

600 

0.71 

11.4 


There remain a number of interesting future directions 
in FACTS device optimization for maximum throughput. It 
will be interesting to include other types of controls, such 
as switching, to determine if the combined control further 
improves system response. It will also be interesting to 
consider problems that are of larger scale than the IEEE 
Polish systems. Finally, future work should also conduct an 
in-depth study of problem structure that results in solutions 
with large gaps between MPF and MFF in order to guide 
the placement of FACTS devices. 
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